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'' The scattering length and effective range of the irS channel are studied in order to 

, characterize the strangeness S = —1 meson-baryon scattering and the yl(1405) resonance. 

QQ . We examine various off-shell dependence of the amplitude in dynamical chiral models to 

' evaluate the threshold quantities with the constraint at the KN threshold. We find that 

the ttS threshold parameters are important to the structure of the 71(1405) resonance and 
provide further constraints on the subthreshold extrapolation of the KN interaction. 

§1. Introduction 

^ . The yl(1405) hyperon resonance has strangeness S = —1 and isospin / = 0, and 

is located just below the threshold of KN. Because yl(1405) is considered to be the 



(N 



quasi-bound state of the KN channel, the structure of 71(1405) is one of the very 
^ . important issues of recent hadron physics, especially to understand ^-nucleon and 

CD ', X-nucleus interactions. At the same time, yl(1405) is a resonance decaying to the 

' ttU channel with / = by the strong interactions. Thus, for the nature of yl(1405), 

ly-^ . dynamics of both ttU and KN is important and gives essential contributions. 

The yl(1405) resonance has been considered as a dynamically generated state in 
meson-baryon scattering for a long time before establishment of QCD, while the pos- 
sibility of the yl(1405) with some three-quark components originated in quark-model 
viewpoint has been also discussed. The first dynamical investigation of yl(1405) in 
terms of meson and baryon degrees of freedom was performed by Dalitz and Tuan.'^' 
^ . There ttU and KN scattering amplitudes were calculated in a coupled-channels 

^ I approach and yl(1405) was obtained as a quasi-bound state of KN. It was also con- 

firmed in a modern point of view that yl(1405) can be predominantly described by 
meson-baryon components in a coupled-channels approach based on chiral dynam- 
ic^ (chiral unitary model), which successfully reproduced yl(1405) in meson-baryon 
dynamics.'^"^ A phenomenological approach also described yl(1405) as a quasi- 
bound state of KN^^ in which an effective interaction of KN was derived so as to 
reproduce the KN scattering length and the mass and width of yl(1405) as 1405 MeV 
and 40 MeV, as suggested by the Particle Data Group (PDG).'Hl' Concerning the 
amplitude above the threshold, this effective potential is essentially consistent with 
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the one extracted from the chiral unitary approachjl^ although the phenomenologi- 
cal model provides quantitatively stronger KN interaction than the chiral potential 
in the region far below the KN threshold. Thus, there are uncertainties about the 
theoretical extrapolation of the KN interaction below the KN threshold. 

Another recent finding is that yl(1405) is composed by two resonance states 
appearing around yl(1405) energies which have different coupling nature to the ttU 
and KN channels .'^''^''SI'I^ According to Ref. Ml, one state which dominantly 
couples to the KN state is located at 1420 MeV instead of nominal 1405 MeV with 
a narrower width, while the other state couples strongly to the ttE channel and 
appears around 1390 MeV with a wider width. Because of the dominant coupling to 
the KN, the higher pole is more relevant to the X-nucleus interaction.'^SI' Thus, it 
is extremely important to reveal the yl(1405) resonance position in the KN channel. 
This double pole structure of yl(1405) also leads to the feature that the yl(1405) 
resonance spectra depend on its production mechanism.!^ '^IS^^ For observation 
of yl(1405) spectra below the KN threshold, an old bubble-chamber experiment in 
K~d — )• 7r~^X!~n gave a hint of the yl(1405) resonance position at 1420 MeV.'22J A 
recent theoretical analysis based on the chiral unitary approachP^ confirmed that the 
71(1405) production is initiated by the KN channel in the K^d — )• vrZ'n, and was able 
to reproduce the yl(1405) production cross section and the shape of the experimental 
observed spectrum. Further detailed experimental information of the ttE spectrum 
will be obtained in ongoing/forthcoming experiments at Jefferson Laboratory,!^ 
GSl^ and J-PARC.E3 

Recent theoretical investigations based on coupled-channels approaches in chiral 
dynamics suggested the double pole structure of yl(1405). Nevertheless, the position 
of the lower pole is dependent on the details of the model parameters, as discussed 
in Ref. I15p . Moreover, Ref. I26p shows that the use of the approximate energy inde- 
pendent potential also changes the pole position drastically. Indeed, since the pole 
position of the higher state is strongly constrained by observed K~p scattering and 
K~ hydrogen experiments, the theories provide very similar pole positions around 
1420 MeV with narrow widths within the framework of several chiral unitary mod- 
els. However, due to lack of ttU — >• ttU experimental data, the properties of the 
lower state are less controlled in theoretical studies. Hence the pole position of the 
lower state has model dependence even in similar models constructed by the same 
concept. The ttE dynamics, including the nature of the lower yl(1405) pole, may 
be important in kaonic nuclear few-body systems. A possible quasibound state of 
KNN was proposed in GO's,!^ and recent theoretical investigations concluded that 
the KNN system has a quasibound state with a large width,'i31'''2Sli'Elli-[32li g^j^^j ^j^g 

theoretical predictions of the energy and width scatter over a wide range. Although 
part of the ambiguity in the model prediction has been attributed to the treatment 
of the three-body dynamic^SU and the energy dependence of the interaction,!^^ it is 
obvious that the present experimental database does not completely constrain the 
KN-ttE interaction. Since the ttU threshold is located at 100 MeV below the KN 
threshold, if K is bound in few-body systems with a large binding energy, the cou- 
pling to the ttU channel is essential to understand the properties of the quasibound 
KNN state. 
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In this paper, we discuss threshold behavior of the ttU scattering^, emphasiz- 
ing its importance in a context of the yl(1405) resonance in chiral dynamics. For 
the physics of yl(1405), it is certainly necessary to have theoretical descriptions of 
ttU and KN dynamics applicable for the energy region between their thresholds, 
which is as wide as 100 MeV. The nature and position of the pole singularity (bound 
state, virtual state, or resonance) around the nU threshold**^! is an important issue. 
Because chiral symmetry indicates that the s-wave ttU diagonal interaction at low 
energy with / = is strongly attractive, a pole singularity may be produced at 
least somewhere, probably close to the threshold, in the complex energy plane. It is 
known that the leading order calculations of the pion-nucleon and pion-pion scatter- 
ing lengths in the chiral perturbation theory work well within 20% accuracy. Thus, 
it is an interesting question whether it is also the case for the scattering length for 
the ttE channel with / = 0. To clarify the nature of the ttU dynamics, we calculate 
the scattering length a-j^x: and the effective range r^ in various models of the ttU 
and KN coupled channels in which different ways of solving Lippmann-Schwinger 
equation are applied in the isospin limit. Finally, we conclude that the ttE scatter- 
ing length and effective range give us new and essential information of the irU-KN 
dynamics between their thresholds and the yl(1405) resonance. We discuss the gen- 
eral relationship between the scattering length and the effective range with the pole 
singularity around the threshold in Appendix. 

§2. Two channel models of KN-ttX! scattering 

We discuss the threshold behavior of the ttU amplitude with / = in connection 
with the yl(1405) resonance appearing below the KN threshold. To describe the nU 
scattering amplitude and the ^1(1405) resonance, it is essential to treat both ttU 
and KN channels simultaneously in a coupled-channels approach. Here we consider 
simple models of coupled KN and ttU channels with 1 = based on chiral effec- 
tive theory. The interaction kernel is given by the Weinberg- Tomozawa interaction 
and we solve coupled-channel Lippmann-Schwinger equation with this interaction 
kernel. In our model, there are two adjustable parameters which come from the reg- 
ularization of the Lippmann-Schwinger equation. These parameters are determined 
by the scattering length with / = 0. We shall also compare different theoretical 
schemes for solving Lippmann-Schwinger equation for the KN-irU coupled system. 

2.1. Framework of the models 

Let us start with Lippmann-Schwinger equation for the t-operator of the ttU 
and KN channels with / = 0: 

i{W) = v + vg{W)i{W) , (2-1) 



*' A preliminary discussion was given in Ref. I33|) . 
**' Here we investigate the coupled- channels scattering amplitudes below the KN threshold. 
Thus, we search the poles in the KN physical sheet of the complex energy plane and classify the 
poles in terms of the Riemann sheets for the nS channel. 
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where W is the scattering energy in the center-of-mass frame, v is the interaction 
kernel, and g is the free two-body Green operator. We solve the Lippmann-Schwinger 
equation in the center-of-mass system. Inserting complete sets of the meson-baryon 
states, we obtain the Lippmann-Schwinger equation for the t-matrix element: 

{h\i{w)\kj) 

where ki and kj are final and initial meson (baryon) momenta with the channel 
indices i, j, respectively. In general, we need off-shell interaction kernel and t-matrix 
element to solve Eq. ()2-2p . In this equation, the two-body meson-baryon state is 
normalized as {ki\kj) = 5ijy^MiN'j{27r)^5^^^ {ki — kj) where Mi = 2uJiEi/Mi with the 



baryon mass Mi, the baryon energy Ei{ki) = \/kf + Mf, the meson energy ijJi{k. 



kf + mf and the meson mass rfij. The on-shell meson (baryon) momentum ki is 
given by 

_ V(iy2 _ (M, + miY){W^ - {Mi - ruiY) 
fe. - ^ • ^^--^^ 

We fix the interaction kernel v based on the chiral effective theory. At the 
leading order of the chiral expansion for the s-wave scattering, we take so-called 
Weinberg- Tomozawa interaction given by 

{h\v\k,) = -^Ui{]ii + %)u, , (2-4) 

where i and j are channel indices of 'kT, {i,j = 1) and KN {i,j = 2), ki is the meson 
momentum for channel i and k^ is fixed to the on-shell energy k^ = uii{ki) in the 
center-of-mass system, is the meson decay constant, Ui is the Dirac spinor for 
the baryon and Cij is the coupling strength which is determined by the flavor SU(3) 
symmetry as 

Note that the meson-energy-momentum dependence of the interaction (j2-4p is the 
consequence of the Nambu-Goldstone theorem, which is manifested as the derivative 
coupling in the effective chiral Lagrangian. Although the flavor SU(3) symmetry 
requires the full coupled channels with r]A and KS, it is found that the effect of these 
channels are small to the scattering amplitude in the energy region of interest.'^ 

With the Weinberg- Tomozawa interaction, which is a contact interaction in the 
coordinate space, we need to regularize the momentum integration appearing in the 
Lippmann-Schwinger equation ()2-2p . As we shall explain in ^2.21 and 12.31 we will use 
two kinds of renormalization schemes. The regularization is done in each channel, 
so that we have two parameters in the models in the form of cut-offs or subtraction 
constants for the irU and KN channels. These parameters should be chosen so as 
to reproduce available experimental data. 
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In this work, we would like to investigate how the ttU threshold behavior is 
constrained by the KN scattering length. For reference, we use a value of the KN 
scattering length with / = 0, ajf^v ~ —1-70 + 0.68i fm, determined by Martin in 
Ref. I34p . This strategy was also taken in constructing "Model (a)" in Ref. [28]). Since 
the scattering length is a complex value, this is enough to determine the two cut-off 
parameters for the ttE and KN channels. This value of the KN scattering length 
leads to the yl(1405) pole position around 1420 MeVW^ 

In the following subsections, we will explain further details of our models, which 
are obtained by different treatments of the regularization of the integral and the 
interaction kernel. These theoretical choices stem from ways of off-shell treatment 
of the scattering amplitudes. In our models, we assume isospin symmetry and use 
isospin averaged masses, 771,^ = 138 MeV, Mx; = 1193 MeV, = 496 MeV and 
Mat = 939 MeV, and the meson decay constant /,r is fixed to be = 92.4 MeV. 
Note that one should utilize the physical masses for a realistic calculation of the 
threshold observables, because the isospin breaking effects may be quantitatively 
large at the threshold.'^S' In addition, for direct comparison to experimental obser- 
vation, Coulomb effects should also be taken into account as done in Refs. [T]). llOp for 
the analyses of KN scattering. However, the present work focuses on a qualitative 
examination of the nature of the ttU channel, and precise evaluations of the ttU 
scattering length and effective range are out of the scope of the present paper. Thus 
we believe that the calculation in the isospin basis using the isospin averaged masses 
without Coulomb corrections is sufficient for the present purpose. 

After obtaining the s-wave scattering amplitudes in our dynamical models, the 
scattering length a and the effective range re in channel i are extracted by 

« =fi{h)\k,^o , (2-6) 

, (2.7, 



dkf \fi{ki) 



where the s-wave scattering amplitude fi{ki) is given by 



m) = -^{kMW)\h) . (2- 



The definition of a and r^ is given with the phase shift 6i by 



1 k"^ 

kiCot6i = - + re^ + -- - . (2-9) 

Qj 



2.2. Model A 



Model A is based on a chiral unitary approach given in Ref. |8|). In this model, 
the interaction kernel for s-wave is given by 



with Cij given in Eq. (j2-5p . According to the N/D method with a simplification of 
neglecting the left-hand cut by = 1, which was developed in Ref.[S]), we can use the 
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on-shell values of the interaction kernel in the Lippmann-Schwinger equation ()2-2p . 
This corresponds to implementing the elastic unitarity. Then the matrix elements 
of V and t inside the g„ and qm integrals in Eq. (|2-2|) can be factorized out from 
the integrals, and one of the momentum integral can be trivially performed because 
the Green function is diagonal. Therefore, we only perform the integral of the free 
Green function. Applying dimensional regularization, we obtain the finite part of 
the Green function as 

d'^qi 2M, 1 



(27r)4 qf -Mf + i€{P- qif - mf + ie 
2Mi ( , , , ^ Mf mf -Mf + W'^ , 



+ ^ [ln(1^2 _ (^2 _ ^2) ^ 2A:iH^) + ln(1^2 ^ (^2 _ ^2^ ^ 

- ln(-T^2 ^ ^j^2 _ ^2) ^ 2kiW) - ln(-VF2 _ (^2 _ ^2^ ^ 2kiW)\ } ,(2-11) 

where = (W^, 0) is the total energy momentum in the CM frame, /i is the scale of 
dimensional regularization, which is to be fixed at = 630 MeV here, and is 
the remaining regularized finite constant term. We call subtraction constant. 

The subtraction constants, d/fiv ^"^^ ^^ttSj are adjustable parameters in this model 
and to be fixed by experimental data or theoretical requirement. Here we determine 
the subtraction constants by the KN scattering length. 

After having done the regularization of the loop integral, we can solve the 
Lippmann-Schwinger equation in an algebraic way: 

2.3. Model B 

Here we explain another potential model (Model B) for the ttU scattering which 
is developed in Ref. [28|) . In this model, the interaction kernel on the mass-shell is 
given by 51 

{h\v\kj) = Fi{ki)X,,{W)F,{kj) = -§|(2VF - Mi - Mj)Fi{h)Fj{kj) , (2-13) 
with 

where Cij is given in Eq. ()2-5p . and Fi{ki) is a separable form factor of dipole type. 
Since we introduce the form factors in Eq (|2-13p . the on-shell kernel has different form 



*-* Since we are using a diilerent normalization of the t- matrix element from Ref. I28p . we have 
some trivial difference in the factor of Eq. (|2-13p than the original form in Ref. I28p. 
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used in Model A. It is, however, noticed that Xij{W) has the same energy dependence 
as the interaction kernel used in Model A except the relativistic correction factor. 

For the off-shell behavior of the interaction kernel, using the same form of the 
separable form factor as in Eq. ()2-13p with off-shell momenta, we define 

{qMW)\q,) = Fi{q,)K,{W)F,{qj) . (2-15) 

The cut-off parameter appearing in each channel will be fixed by the inputs later. 
We call the model with this energy dependent interaction by "Model B E-dep". 

We also consider an energy independent interaction given by taking the 0th 
components of /cf and kj of Eq. ()2-4p and the static approximation: 

C- ■ 

{ki\v\kj) = -—^{u}i{ki) +ujj{kj))Fi{ki)Fj{kj) 

J IT 

- -77l("^* + mj)Fi{ki)Fj{kj) . (2-16) 

These kinds of the energy independent potential are often used in few-body cal- 
culations and are favorable for non-relativistic quantum calculations. The off-shell 
behavior is defined in the same way as Eq. ()2T5p . We call the model with this 
energy-independent interaction by "Model B E-indep". This energy-independent 
potential may contradict the energy expansion scheme in chiral effective theory at 
energies far from the threshold. Nevertheless, it is interesting to compare the results 
with this potential and others, in order to see the importance of energy dependence 
of the potential. 

Since we introduce the form factor in the interaction kernel, the integration is 
regularized: 

^ / /^'^f3!^ {k,\v\q„,){q^^^^ 



n=l 



with 

{h\i{W)\kj) = F,{h)nj{W)Fj{k,) , 

with En = \/qn~\~ ^'^'^ = ■\/'ln~^ '^n- (|2T7p , we use the fact that the t- 

matrix elements can be factorized in separable forms when the separable interactions 
are taken into account. The i-matrix element is obtained again in an algebraic way 

as 

where 

Gn{W) ^ [ FMFM ^2-19) 



8 



Y. Ikeda, T. Hyodo, D. Jido, H. Kamano, T. Sato and K. Yazaki 



§3. Threshold parameters in ttX! channel and A(1405) 

In this section, we show the results of the ttE scattering length and effective range 
obtained by our two-channel models. First of all, we fix the parameters appearing in 
each model so as to reproduce the KN scattering length a/f^ = — 1.70 + 0. 68i fm.'SS 
To characterize the subthreshold behavior of the amplitude, we also calculate the 
pole positions of the yl(1405) resonance in the complex energy plane. The poles are 
searched in the KN physical Riemann sheet and classified in terms of the two ttU 
sheets; bound states in the first Riemann sheet, while virtual states and resonances 
are in the second Riemann sheet. 

Let us first show the result of Model A with dimensional regularization of the 
loop function. We find two solutions to reproduce the KN scattering length (a^^v ~ 
-1.70 + 0.68i fm); we shah cah these two solutions "Model Al" and "A2": 

dnE = -1-67, d^N = -1-79, (Model Al), (3-1) 

dnE = -2.85, d^N = -2-05, (Model A2). (3-2) 

Note that the natural value of the subtraction constant^ is dj^j^ = —1.95 and dj^x: = 
— 1.60 at this scale, so that Model Al is more natural from the theoretical point of 
view. 

With the subtraction constant p-ip . the nU scattering length and effective range 
are obtained as 

a^i: = 0.934 fm, re = -5.02 fm, (Model Al). 

The pole singularities in the scattering amplitudes are found at 

z = 1422 - 16i MeV, z = 1375 - 72i MeV, (Model Al). 

These poles correspond to resonance states located between KN and ttE thresholds. 
The pole positions obtained from Model Al are similar to the ones obtained in 
standard parameterizations of chiral unitary models; the KN bound state and the 
ttS resonance .f^'fl^ 

The other solution (Model A2) in Eq. ()3-2p provides the ttU scattering length 
and effective range as 

a,ri: = -2.30 fm, rg = -5.89 fm, (Model A2), 

and we find the pole singularities at 

z = 1425 - 111 MeV, z = 1321 MeV (bound state), (Model A2). 

In this case, one pole appears at z = 1425 — Hi MeV which may be interpreted as 
the yl(1405) resonance, while the other pole is found in the first Riemann sheet of the 
nU channel below the threshold, corresponding to a bound state of ttE. Following 
the discussion in Ref. |2|), the large negative value of the subtraction constant is 
equivalent to the enhancement of the interaction strength. Thus, this parameter set 
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provides the ttU interaction stronger than Model Al and the scattering length in the 
irE channel becomes negative. A general discussion of the relation of the scattering 
length and pole position of scattering amplitude is given in Appendix. 

In Model B, we utilizes the separable potential with form factors and have two 
options of the energy dependence of the interaction kernel ()2-13p and (j2-16p . The 
KN scattering length a/f^ = — 1.70 + z0.68 fm is reproduced by the following cut-off 
parameters: 

A^s = 1005 MeV, yl^^ = 1188 MeV, (Model B E-dep), (3-3) 
A^s = 1465 MeV, A^j^ = 1089 MeV, (Model B E-indep). (3-4) 

With these parameters, the ttE scattering length and effective range are obtained as 

= 1.44 fm, = -3.96 fm, (Model B E-dep), (3-5) 

= 5.50 fm, rg = -0.458 fm, (Model B E-indep). (3-6) 

and the pole positions are 

z = 1422 - 22i MeV, z = 1349 - 54z MeV, (Model B E-dep), (3-7) 

z = 1423 - 29i MeV, z = 1325 MeV (virtual state), (Model B E-indep). (3-8) 

The result of Model B E-dep is very similar to that of Model Al, although different 
regularization schemes are applied. For the energy independent interaction case, 
one of the poles is obtained as a virtual state located in the second Riemann sheet 
of the ttU channel below threshold. Thus, the scattering length is to be a large 
positive value. This is also understood by the comparison of the cutoff parameters 
in the model. Since the large cutoff parameter effectively enhances the interaction 
strength, the nU attraction in Model B E-indep is stronger than that of Model B 
E-dep, leading to the formation of a virtual state below the threshold. It should be 
noted that Model B E-indep also develops two poles as shown in Ref. [26j) . although 
one of them is a virtual state. 

We would like to emphasize that all the parametrizations are constrained by the 
KN scattering length very well, but ttE scattering lengths and pole positions in these 
parametrizations are quite different. This means that only with the KN threshold 
behavior, the nS and KN coupled-channels amplitudes cannot be determined, even 
though the interaction kernels are fixed by theoretical consideration. This is clearly 
seen by plotting the scattering amplitude / as functions of W (Fig. [1]). In the left 
panel, the KN amplitudes around threshold coincide with each other since we use 
the KN scattering length to fix the model parameters. The agreement of the two 
models holds down to ~ 1420 MeV, but below that energy, the extrapolation 
cannot be controlled by the KN scattering length. The difference is significant in 
the right panel where the amplitude of the ttE channel is plotted. Of course, we 
have the experimental data of the ttE mass spectrum and the threshold branching 
ratios K~p — ?• ttU and ttA observed in decay of kaonic hydrogen,'3Sl''[3Zli which should 
reduce the difference between models. It is nevertheless important to keep in mind 
that the KN scattering length does not fully constrain the structure of the yl(1405) 
resonance and the KN amplitude far below the threshold. 
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1320 1360 1400 1440 1320 1360 1400 1440 

W [MeV] W [MeV] 

Fig. 1. Scattering amplitudes in KN (left panel) and nS (right panel) channels. The thick (thin) 
curves represent the results from model Al (A2), and real and imaginary parts are plotted by 
solid and dashed curves in both panels. 

§4. Discussion 

We have investigated the scattering length and effective range of the ttE channel 
with / = using the two-channel models of the ttU and KN coupled systems. We 
have considered two kinds of models with different ways to solve the Lippmann- 
Schwinger equation. They are constrained by the KN scattering length, in order 
to see how these inputs constrain the ttE threshold behavior and the yl(1405) pole 
positions. The results are summarized in Table [H 

First of all, we find that the KN scattering length constrains the position of 
the pole near the KN threshold, which we call the (higher) yl(1405) pole, well 
around 1420 — 20i MeV in all of our models. Thus, the value of the KN scattering 
length with / = can be one of the important quantities to fix the pole position 
of yl(1405) which strongly couples to the KN channel. On the other hand, the ttU 
scattering length and effective range are obtained with very different values. This 
is a consequence of different predictions of the lower pole positions. This means 
that the KN scattering length alone cannot constrain the scattering amplitude at 
far below threshold. In contrast, the ttE scattering length and effective range are 
sensitive to the lower pole position, which will give important constraints on the 

Table 1. Summary of numerical results. The model parameters are determined so as to reproduce 
a value of the KN scattering length with 7 = 0, a^jv = —1.70 + 0.68i fm. In each model, pole 
1 is found as a resonance located between the ttS and A'A*" threshold. For pole 2, (R), (B) and 
(V) denote resonance, bound state and virtual state in the nS channel, respectively. 



Model 


Al 


A2 


B E-dep 


B E-indep 


parameter (nS) 


d-^E = —1.67 


d^s = -2.85 


A^s = 1005 MeV 


A^s = 1465 MeV 


parameter {KN) 


dji^ = -1.79 


rfjj^r = -2.05 


yljfiv = 1188 MeV 


Aiiff = 1086 MeV 


pole 1 [MeV] 


1422 - 16i 


1425 - lli 


1422 - 22i 


1423 - 29i 


pole 2 [MeV] 


1375 - 72i (R) 


1321 (B) 


1349 - 54i (R) 


1325 (V) 


a^s [fm] 


0.934 


-2.30 


1.44 


5.50 


Te [fm] 


-5.02 


-5.89 


-3.96 


-0.458 


aj^pf [fm] (input) 


-1.70 + 0.68i 


-1.70 + 0.68i 


-1.70 + 0.68i 


-1.70 + 0.68i 
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Table 11. Summary of numerical results in the tvU single calculations by switching off the KN 
channel effects. The model parameters are same as Table U 



Model 


Al 


A2 


B E-dep 


B E-indep 


pole 2 [MeV] 


1366 - 88i (R) 


1331 (V) 


1345 - 79i (R) 


1268 (V) 


a^s [fm] 


0.77 


34.4 


0.935 


1.42 


r-e [fm] 


-4.59 


-4.89 


-3.58 


-0.831 



subthreshold extrapolation of the KN amplitude. In models Al and B E-dep we 
use essentially the same interaction kernel but take different off-shell dependence. 
The results are qualitatively the same but there are quantitative differences in the 
irU threshold parameters. Thus, we find that the additional ttE threshold quantities 
can be useful also to constrain dynamics of the KN-ttU coupled channels below the 
KN threshold. 

In the model A2 and the model B with energy independent potential, the ttU 
interaction is so strong that they give a bound state and a virtual state, respectively. 
In both cases, the ttU attraction is effectively enhanced by the cutoff parameter. 
Though the presence of such a bound or virtual state is in contradiction with the 
result of the more refined calculations with the chiral unitary approaches, which 
take account of the K~p scattering data, it is important to clarify the position of 
the lower yl(1405) pole by experimental observation. This can be done by observing 
the sign and order of magnitude of the ttU scattering length. If the ttU scattering 
length in the / = channel would have a negative value, there could be a bound 
state of TT and U with / = 0. If the scattering length is positive with as large a value 
as 5 fm, there is a virtual state below and close to the ttU threshold. This also shows 
the relevance of the ttU threshold behavior for the subthreshold extrapolation of the 
KN amplitude. 

It is also interesting to mention that the present model provides the higher 
yl(1405) pole around 1420 — 20i MeV irrespectively of the off-shell behavior and the 
subthreshold extrapolation to much lower energy region, when we use the jVIartin's 
value of the KN scattering length with I = to fix two parameters appearing 
in the unitarization procedure. This fact was implicitly shown in Ref. [28|) . Since, 
in our simple two-channel model, the interaction is fixed by chiral effective theory, 
we would say that this is a consequence of the unitarization of the leading order 
chiral interaction. Nevertheless, for the definite conclusion on the pole position, it 
is certainly necessary to refine the analysis of IVIartin done in 80's by including new 
measurements of K~p, such as KEK, DEAR and forthcoming SIDDHARTA data, 
and also to use more sophisticated models, for instance, including other channels, 
especially 1 = 1 amplitudes and isospin breaking effects. 

In order to see the ttE-KN coupled-channels effects, we calculate the scattering 
length, the effective range and the pole position in the ttU single channel by switch- 
ing off the KN channel. We use the same model parameters as before, which are 
determined by the KN scattering length in the full coupled channel. The results are 
summarized in Table [III which shows that the coupled-channels effects contribute 
moderately as an attractive interaction to the ttE channel at the threshold, which 
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Table III. Summary of numerical results with the /l(1405) pole position being 1406.5 — 25i MeV. 
In each model, pole 1 is found as a resonance located between the irS and KN threshold. For 
pole 2, (R), (B) and (V) denote resonance, bound state and virtual state in the nS channel, 
respectively. 



Model 


Al' 


A2' 


B' E-dep 


B' E-indep 


parameter (nS) 


d,rs = —1.66 


d^s = -2.68 


A^s = 860 MeV 


A^E ^ 1070 MeV 


parameter {KN) 


dj^jv = —1.93 


dn^ = -2.44 


yljfiv = 1290 MeV 


^KN = 1175 MeV 


pole 1 [MeV] (input) 


1406 - 25i 


1407 - 25i 


1405 - i25 


1406 - i26 


pole 2 [MeV] 


1381 - 65i (R) 


1304 (B) 


1351 - i64 (R) 


1306 (V) 


a^s [fm] 


1.03 


-0.75 


1.38 


2.77 


re [fm] 


-5.23 


-14.32 


-3.90 


-0.509 


aRN [fm] 


-1.37 + 0.4H 


-0.96 + 0.24i 


-1.48 + i0.38 


-1.67 + i0.40 



can be seen in the values of the scattering length. The scattering lengths in the 
model Al and B E-dep slightly change, while in the model A2 and B E-indep the 
coupled-channels effects are largely seen in the scattering length. In any cases, the 
pole position is consistent with the classification of the poles in terms of (a, rg) given 
in Appendix. 

As mentioned above, the scattering length a/^jy = —1.70 + zO.68 fm leads to 
the pole of yl(1405) around 1420 MeV. Although the pole position is not a direct 
observable, the obtained value deviates from the nominal value shown in PDG,'^^!' 
1406.5 — 25i MeV. To examine a different input, we set up the models so as to 
reproduce the pole position of the PDG value, and perform the same analysis. In 
this case, we again find two solutions in model A (Al' and A2') and model B (B' 
E-dep and B' E-indep). Thus, we have four solutions as summarized in Table Hill 
which provide different values of the scattering length of the KN channel. This 
implies that, if yl(1405) is located at as deeply as 30 MeV below the KN threshold, 
the yl(1405) pole position is also affected by the ttU dynamics and the pole position 
alone cannot constraint the KN threshold quantities. Again, precise information of 
the ttU dynamics is necessary to understand the structure of 71(1405) further and 
to have models constrained more. It is also important mentioning that, even in 
such two-body ttU and KN dynamics, both of these channels are important and 
a quasibound state with 30 MeV binding energy is generated by a consequence of 
both ttU and KN dynamics. For investigation of more deeply bound systems, which 
would be seen in KNN, since the irE dynamics should become more important in 
such systems, one needs experimental information of the irU interaction and more 
serious consideration of the irE channels. 

It is also interesting to compare the present results with those obtained in more 
refined models. In table IIV|, we show the ttU scattering length and effective range 
calculated in various chiral unitary approaches together with the pole positions of 
yl(1405). We adopt the models with the Weinberg- Tomozawa interaction in Refs.[7I), 
[5|).ll0 p .[TT |) where all the channels consisting of the octet baryons and the octet 
pseudoscalar mesons with S = —1 and / = {ttU, KN, r]A and KS) are included 
and the model parameters are determined so as to reproduce the total cross sections 
of the K~p scattering and the threshold branching ratios of K~p — )■ Ytt with Y = A 
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Table IV. Scattering lengths and effective ranges of the nS channel with 7 = in chiral unitary 
approaches and the phenomenological potential model (EAY). The value of the pole positions 
of yl(1405) and the result of KN scattering length with J = in each model are also listed. The 
calculation is performed with isospin symmetric masses. 



model 






BN-#i2J 






pole 1 [MeV] 


1427 - 17i 


1428 - 17i 


1434 - 18i 


1421 - 20i 


1409 - 20i 


pole 2 [MeV] 


1389 - 64i 


1400 - 76i 


1388 - 49i 


1440 - 76i 


1299 (V) 


a-^E [fm] 


0.789 


0.693 


0.770 


0.517 


2.24 


Te [fm] 


-6.27 


-6.50 


-7.18 


-8.63 


-0.447 


a-KN [fm] 


-1.65 + 0.89i 


-1.63 + l.Oli 


-1.49 + l.lli 


-1.45 + 0.872 


-1.76 + 0.42i 


U [MeV] 


103.7652 


106.95 


111.2 


120.9 





or E in the observed strong decay of kaonic hydrogen. For comparison with our two- 
channel model, as in Ref.llSp. we use the isospin averaged masses with the subtraction 
constants and the pion decay constant being kept fixed to be the original values in 
these references. These results are very similar to our models Al and B E-dep. The 
results obtained by the chiral unitary models are qualitatively consistent with each 
other. Nevertheless, in detail, the positions of the lower pole of yl(1405) and the 
ttE scattering lengths are differently predicted. When the pole is closer to the ttU 
threshold (ORB, BNW), the nX! scattering length is relatively large. It is also worth 
mentioning that the parameters of BIVIN were determined by a fitting using the 
currently available KN scattering data and Ref. ITT]) found a large uncertainty of the 
position of pole 2 in their fitting. Thus, the quantitative determination of the ttU 
scattering observables could be a good guidance to refine the models for the scattering 
amplitude further and to locate the position of the lower pole of 71(1405). We also 
show, in Table llVt the result obtained by the phenomenological potential modeO 
developed in Ref.lSSp. in which irXJ-KN coupled channels with / = were considered 
with an energy-independent potential. The interaction range and strengths of KN- 
KN and ttE-KN were determined by K~p scattering data, kaonic atom data and 
the pole position of yl(1405) at 1405 — 20i IVIeV,'^ while the interaction strength of 
irU-irX! was fixed so that the ratio of the whole potential strengths of irX! and KN is 
to be 4/3 as motivated by the chiral interaction, which is, however, slightly different 
from the interaction strength in the chiral model given in Eq. (j'2-16p . The factor 
Cij in Eq. ()2-16p does have the ttE/KN ratio equal to 4/3 according to Eq. (|2-5p . 
(These two potential strengths coincide in the flavor SU(3) limit.) The result of the 
phenomenological model is found to be very similar to that obtained in B' E-indep 
given in Table lllli This is because the phenomenological potential model uses the 
energy-independent interaction and the parameters were fixed by the yl(1405) pole 
position at 1405 — 20i IvleV. In addition, the ttU interaction strength is similar to 
the chiral models. These are also the case in B' E-indep. A virtual state is found 
also in the phenomenological potential model due to the attractive interaction in the 

*■* In this calculation we have used the isospin averaged masses given in page[5]and obtained the 
pole position at 1409 — 20i MeV as shown in Table [TVl It seems that in Ref. I38[l the value of the KN 
threshold is given by sum of the proton and K~ masses. With these values, we have reproduced 
the yl(1405) pole position at 1405 - 20i MeV. 
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ttE channel. 

The ttU channel is the lowest energy meson-baryon state with / = and 5 = 
— 1, and it is believed that there is no bound state below the threshold. One may 
consider that such a lowest state with one pion may be well described by the chiral 
perturbation theory, as in the case of the irN system. Using the Weinberg- Tomozawa 
interaction (|2-4p with Ct^s = 4, we obtain the scattering length and effective range 
at the leading order of the chiral perturbation theory as 

a^E = = 0A55 fm, (4-1) 

..._^/>i±l?£|)._4.52ta. ,4.2, 

with the -kE reduced mass /ijri; = fn-,^M^/{mT^ + M^:)- Comparing the results ob- 
tained by the chiral perturbation theory with those calculated in dynamical models, 
we find that the deviation of the results of the dynamical models from those of the 
chiral perturbation theory is much larger than 20%, in contrast to the vrA^ scatter- 
ing case where the perturbative calculation works well. This is understood by the 
existence of the pole singularity around the threshold. The interaction strength in 
the 7ri7 channel is C^^s = 4, which is close to the critical coupling constant to have 
a bound state.'^SI' The strong attraction generates a pole around the threshold and 
hence the perturbative calculation may not be applicable. On the other hand, the 
coupling constant in the vrA^ channel is CttAt = 2 {Ct^m = —1) for / = 1/2 (/ = 3/2), 
which is weaker than the vrZ" case (repulsive) so that no pole singularity appears 
around the threshold. In this way, the determination of the vrZ" threshold quantities 
with 1 = gives us an interesting test for the applicability of the chiral perturbation 
theory. 



§5. Conclusion 



We have investigated the vrZ" scattering length and effective range in various 
models, solving Lippmann-Schwinger equation for the ttE and KN coupled channels 
with the interaction kernel obtained by the leading order chiral perturbation theory. 
We have used different regularization schemes to integrate the Lippmann-Schwinger 
equation, which is equivalent to applying different treatments of off-shell scattering 
amplitudes. We find that the constraints from the KN scattering length determine 
the theoretical descriptions of the scattering amplitude around KN threshold and the 
higher pole position of yl(1405), whereas it is not enough to fix the precise position 
of the lower energy pole and the scattering amplitude around the ttE threshold. 
Consequently, the vrZ" scattering length and the effective range will give us new 
additional information for the ttU-KN dynamics and the structure of yl(1405). 

At this moment, there is no experimental information of the ttE threshold quan- 
tities. Although scattering experiments are impossible, we may utilize the interfer- 
ence of the final state interaction in the decay of the heavy particle,'3Sl' in analogy 
with the determination of the vrvr scattering length in Ref. HT]) . Moreover, recent 
development of lattice QCD may bring the information of the scattering length.SS' 
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As a first step, the sign of tlie scattering length can be used to check the existence 
of the bound state below the threshold. In the long run, the determination of the 
magnitude of the scattering length and the effective range will provide the detailed 
information of the pole structure of 71(1405) and the severe constraints for the ttE- 
KN scattering amplitude. We hope that the present semi-qualitative analysis on 
the ttU scattering length and effective range gives a good guideline for forthcoming 
experiments and lattice QCD calculations. 
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Appendix A 

Scattering length and effective range 

Let us consider an s-wave single-channel two-body scattering amplitude f{k) 
with momentum k. This is given in terms of the S matrix element with / = by 

^^''>~ 2ki ~kcot6-ki' 

where 6 is the s-wave phase shift. The scattering length a and the effective range r^ 
are defined as the expansion coefficients of k cot 5 in terms of k'^ around k = 0\*^\ 

kcot6 = - + re— -\ . (A-2) 

a 2 

Taking only the first two terms in Eq. (|A-2p . we write the scattering amplitude as 

This amplitude has a pole, when the momentum k satisfies the condition 

jk"^ -ki + ^=0 . (A-3) 

*-* The sign of a is opposite to the standard convention in nuclear physics. In our convention, 
f{k) — !> a when fc 0. 
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If Tg = 0, the solution of Eq. ()A-3p is given by 

k = -- . (A-4) 

a 

If there are no open channels below the threshold {k = 0), the scattering amplitude 
is real at A; = 0, and consequently the scattering length a is also real. In this case, 
the pole singularity given by Eq. ()A-4p lies on the imaginary axis on the complex k 
plane. The solution with Im A; > (Im A; < 0) corresponds to the bound (virtual) 
state, so we have a bound state for a < and a virtual state for a > 0. When the 
bound state is formed, the scattering length looks as if the interaction is repulsive. 

Next, we study the case of ^ 0. In this situation, the inverse of the scattering 
amplitude is given by a quadratic function, so that we always have a pair of poles 
for given a and except for the case of two- fold root. The solutions of Eq. ()A-3p 
are 



k=—±—\ (A-5) 

Te Te V a 

Defining — 2re/a — 1 = D, the solutions can be classified by the sign of D: 

• When D < 0, the solutions are pure imaginary. 

• When D > 0, the solutions are complex with the same imaginary part. 

The case D < corresponds to a bound state or a virtual state without width, 
while a virtual state with finite width or a resonance can be formed for the case 
D > 0. For later convenience, we schematically illustrate the pole position and the 
resonance/virtual/bound state in Fig. [5J The poles with Im A; > (Im A; < 0) are 
mapped onto the first (second) Riemann sheet of the energy plane. As mentioned, 
there are always two poles for given (a, rg), but the scattering amplitude above the 
threshold is mainly affected by the pole close to the scattering region Re A; > 
(triangles in Fig. [2]). In the following, we classify the solutions by the character of 
the pole most relevant to the amplitude, which appears at larger Im k for the state 
without the width or at Re A: > for the state with the width. 

Let us first consider the D < case. The parameter region for D < is given 

by 

fa>0, rg>-a/2 
[a < 0, re < -a/2 

The pole positions are 

Re A;poie =0, Im A;poic = ^ (^+1 ± 

In this case, both the two solutions lie on the imaginary k axis. We can further 
classify the solutions by the values of re and D: 

re > 0, \/—D > 1, virtual state without width 

re > 0, \/—D < 1, bound state , . , 

(A-7) 



< 0, \J —D > 1, bound state 

< 0, \/—D < 1, virtual state without width 
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(bound state) 



Re 



(resonance state) 



X 



X 



w {virtual state) 



Fig. 2. Schematic illustration of the pole position and the resonance/virtual/bound state in the 
complex k plane. Triangles represent the poles which are relevant to the amplitude above the 
threshold, and the crosses denote the other poles. The dotted line (Re -Epoic ~ 0) is the boundary 
between the regions of resonances and virtual states. 



Note that > 1 means 2re/a > 0, which corresponds to the case where a and 

Te has the opposite sign. Thus, in terms of a and re, we have 



{a > 0, Te > —a/2, virtual state without width 
a < 0, ?*£ < —a/2, bound state 

Next we consider the D > case. The parameter region is given by 



(A-8) 



(A-9) 



ia> 0, < -a/2 
a <0, Tg > -a/2 

The real and imaginary parts of the pole can be extracted from Eq. ()A-5P as 

1 1 
Re /cpoie = ± —Vd, Im kpoic = —■ 
re ' re 

In this case, the pair of poles should emerge in the symmetric position with respect 
to the imaginary k axis. If rg > 0, two poles with finite width appear in the first 
Riemann sheet of the complex energy plane. Since such poles are forbidden by 
the causality, we consider them to be an artifact of the truncated effective range 
expansion and expect them to disappear when higher order terms in k are included. 
We therefore classify the region of {a, re) giving rise to such poles as "no state". 
From Eq. ()A-9|) . we find that this is the case for a < 0. For rg < 0, the two poles 
are interpreted as a virtual state with width or a resonance state, depending on the 
ratio of the real and imaginary parts. The boundary in the complex k plane for 
the virtual and resonance state (with nonrelativistic kinematics -Bpoic = ^poie/^^) 
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given by Re -Epoio = in the second Riemann sheet, namely 

{Im /cpoic < 0, I Re /cpoic|/|Ini A;poic| > 1; resonance state 
Im /cpoic < 0, |Re /cpoic|/|Im /cpoid < 1, virtual state with width 

The classification of the solution is then given by 



(A-IO) 



D > 1, resonance state 

D < 1, virtual state with width 



(A-ll) 



Thus, combining (|A-9p and (|ATip . the parameter space for obtaining resonance/virtual 
state is given by 



a > 0, Tg < — a, resonance state 

a > 0, —a/2 > > —a, virtual state with width 

a < 0, Te > —a/2, no state 



(A-12) 



We summarize Eq. ()A-8p and ()A-12p in a-r^ plot in Fig. [31 If a is negative, there is 
a bound state. For positive a, the ratio of a and rg determines the location of the 
singularities. Setting = 0, we recover the classification by the scattering length 
Eq. (|X4l) . 



1 

(no state) 


(virtual state ■ 
without width) I 


(bound state) 


^rn^ ► 

4 

\ (virtual state 
'\ with width) 

(resonance state) \ ^ 



?• = - a/2 



Fig. 3. Parameter region corresponds to resonance/virtual/bound state in a-r^ plot. 

To illustrate the relation between the threshold parameters and pole position of 
the amplitude, we solve the single-channel version of the model shown in section [221 
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d [dimensionless] Re Z [MeV] 

Fig. 4. The nS threshold quantities (left) and the pole position of the scattering amplitude (right) 
with various subtraction constants d. All symbols are plotted at each 0.5 step of the subtraction 
constant. In the left panel, the effective range is plotted with opposite sign and a/2 is also 
shown for positive scattering length, for convenience of the comparison with Eq. (|A-12|) . In 
the right panel, there are always two poles in the complex energy plane, but we only show the 
most relevant pole to the observables defined on the real axis above threshold. The poles on the 
first (second) Riemann sheet are plotted by the crosses (triangles), which are slightly shifted to 
positive (negative) direction of the imaginary axis when they are on the real axis. 

We choose the ttU channel where we obtain a resonance state with the subtraction 
constant d = —2. The scattering length and effective range are 

a = 1.06 fm, Te = —4.67 fm, 

in accordance with the resonance solution in (|A-12p . r^ < —a < 0. When we increase 
the absolute value of the subtraction constant, the interaction strength is effectively 
enhanced,'^' so the pole is expected to become a virtual state and eventually a bound 
state. 

We calculate the scattering length, the effective range and the pole positions 
by varying the subtraction constant d from —2 to —3.5 as shown in Fig. HI The 
effective range is shown with opposite sign for comparison. As the absolute value of 
the subtraction constant is increased, both the scattering length and effective range 
(—re) become large. The resonance pole moves to the lower energy region with 
decreasing the width. This indicates the effective enhancement of the interaction 
strength. 

In contrast to the mild d dependence of the effective range, the scattering length 
increase rapidly for d < —2.6. At d = —2.7, the scattering length becomes larger 
than the effective range with negative sign. At the same time, the resonance pole 
moves below the threshold and becomes a virtual state with finite width. When 
the subtraction constant reaches d = —2.8, the half of the scattering length exceeds 
the effective range a/2 > — rg. This coincides with the appearance of the virtual 
state pole on the real axis. As we further increase the magnitude of the subtraction 
constant, the pole finally becomes a bound state for d < —2.9. Before that, the scat- 
tering length becomes very large positive value. Once the bound state is formed, the 
scattering length is obtained as a negative value whose absolute value is decreasing. 
These behaviors are fairly consistent with the classification given in Eqs. (|A-8p and 
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